Measurement of Analyzing Power for Proton-Carbon Elastic Scattering in the 
Coulomb-Nuclear Interference Region with a 22-GeV/c Polarized Proton Beam 
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The analyzing power for proton-carbon elastic scattering in the coulomb-nuclear interference re- 
gion of momentum transfer, 9.0 x 10 -3 < — t < 4.1 x 10 -2 (GeV/c) 2 , was measured with a 21.7 GeV/c 
polarized proton beam at the Alternating Gradient Synchrotron of Brookhaven National Labora- 
tory. The ratio of hadronic spin- flip to non-flip amplitude, rg, was obtained from the analyzing 
power to be Rer 5 = 0.088 ± 0.058 and Imr 5 = -0.161 ± 0.226. 

PACS numbers: 25.40.Ve, 13.88.+e, 13.85.Dz, 29.27.Hj 



The analyzing power, An, defined by the left-right 
asymmetry of the cross sections in the scattering plane 
normal to the beam polarization, can only arise from the 
interference between a spin-flip and a non-flip amplitude 
and thus provides important information on the spin- 
dependence of the interactions. In high energy proton- 
proton (pp) jl], D and proton-nucleus (pA) || |4j] elastic 
scattering at very small momentum transfer, An origi- 
nates from the interference between the electromagnetic 
(coulomb) spin-flip and the hadronic (nuclear) non-flip 
amplitude. This analyzing power has been calculated 
to have a maximum value of about 4 % around mo- 
mentum transfer — t = 2 x 10~ 3 (GeV/c) 2 , and to de- 
crease as \t\ increases. An in this coulomb-nuclear in- 
terference (CNI) region of momentum transfer has been 
suggested as a sensitive probe of the hadronic spin-flip 
amplitude g. The hadronic spin- flip amplitude is char- 
acterized by the ratio of the spin-flip to the non-flip am- 
plitude, r 5 — (rn/y/^i) x F^/lmF^, where F^ s are the 
hadronic parts of the spin- flip amplitude, F s , and non- 
flip amplitude, F , and m is the nucleon mass. Exis- 
tence of the hadronic spin-flip amplitude introduces a 
deviation from An calculated with no hadronic spin-flip 
amplitude. There are theoretical estimates and experi- 
mental evidence for the possibility of a non-zero hadronic 
spin-flip amplitude reflecting helicity non-conservation in 
the direct channel of the reaction Q . A hadronic spin-flip 
amplitude non-vanishing at high energies carries impor- 
tant physics information on static properties and on the 
constituent quark structure of the nucleon. In the frame- 



work of Regge phenomenology, hadronic amplitudes sur- 
viving in the high energy asymptotic region are described 
by Pomeron exchange. The Regge pole model does not 
limit the value of the Pomeron spin-flip amplitude, but 
leads to a vanishing contribution to An due to equal 
phases of the spin-flip and non-flip parts of the Pomeron 
exchange. An in the CNI region remains sensitive to the 
Pomeron spin-flip amplitude because of the phase dif- 
ference between the hadronic and coulomb amplitudes, 
which is close to tv/2. As was pointed out by Kope- 
liovich and Trueman j|, the usage of a carbon target 
(an isoscalar nuclear target) has an important advantage 
of eliminating the contribution of the isovector Reggeons 
and thus allows one to probe the Pomeron spin-flip am- 
plitude through An at medium high energies. 

Measurement of An for proton-carbon (pC) elastic 
scattering in the CNI region has been proposed not only 
to investigate a hadronic spin-flip amplitude, but also 
to provide a method for high energy proton polarime- 
try at the Relativistic Heavy Ion Collider (RHIC) of 
Brookhaven National Laboratory (BNL) Predicted 
properties of An (large cross section and the weak y/s- 
dependence), make this process ideal for beam polariza- 
tion measurements at RHIC. 

The E950 experiment at the Alternating Gradient Syn- 
chrotron (AGS) of BNL was carried out to measure An 
for pC elastic scattering in the CNI region at 21.7 GeV/c. 
The measurement was made in the AGS ring using an in- 
ternal carbon target. We identified pC elastic scattering 
by detecting only recoil carbon ions, which had energies 
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FIG. 1: Schematic layout of the E950 experimental setup (not 
to scale). 



ranging from 300 keV to 2000 keV and recoil angles of 
89 — 90° with respect to the beam direction. The for- 
ward scattered proton, at very small scattering angles, 
was not detected. In this letter, we present the result of 
the measurement of An and the constraint on r$ . 

The transversely polarized proton beam was acceler- 
ated in the AGS with the beam polarization preserved by 
a partial Siberian snake j?J and an rf dipole || . A fiat- 
top of 520 ms in the AGS accelerating cycle was used for 
the experiment. The polarization direction was flipped 
for each AGS cycle to reduce systematic uncertainties in 
the An measurement. The beam was bunched so that 
we were able to determine the timing of a beam crossing 
at the carbon target and separate slow recoil carbon ions 
from prompt backgrounds. The AGS beam consisted of 
one bunch of about 6 x 10 9 protons with a width of ~25 ns 
and with a rotation period of 2.7 /is. 

The beam polarization was measured concurrently by 
experiment BNL-AGS E925, based on the analyzing 
power of pp elastic scattering at 21.7 GeV/c and at 
-t = 0.15 ± 0.05 (GeV/c) 2 , A% = 0.0400 ± 0.0048 §. 
For each AGS cycle, the beam was debunched and was 
extracted to the CH2 target of the polarimcter after the 
data taking time window for the E950 experiment, pp 
elastic scattering was identified by detecting both scat- 
tered and recoil protons. The average beam polariza- 
tion during the E950 experiment was determined to be 
V B = 0.407 ± 0.036(stat) ± 0.049(syst). 

Carbon micro-ribbon targets flC| ] were employed by 
bundling several ribbons together (each of which were 
6 /im wide x 3 cm long by 3.7±0.2 /ig/cm 2 ) in order 
to get sufficient counting rate. The effective target thick- 
ness, needed for carbon recoil energy loss corrections, was 
estimated to be 6.6±4.0 /xg/cm 2 from the multiple scat- 
tering width of the recoil carbon angular distributions. 

The detector system consisted of two symmetric left 
and right arms detecting recoil carbon ions. Fig. [TJ shows 
a schematic layout of the experimental setup. Each de- 
tector arm consisted of a thin 5 /ig/cm 2 carbon foil as a 
secondary electron emitter, an accelerating grid, an elec- 
trostatic mirror and a micro-channel plate (MCP) assem- 
bly, all followed by a silicon strip detector (SSD). The 
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FIG. 2: (a) The correlation between the energy deposit, E, 
and the timing, Tssd, in one of the SSD strips. Recoil carbon 
events are seen in the distribution at longer times. The flat 
band at shorter times is from background, which is also seen in 
target-empty data, (b) The mass distribution in the measured 
i-range. Events selected with the 2.5-sigma mass resolution 
cut are shown as a hatched area. 



SSD had a thickness of 400 /im, sensitive area of 24 mm 
wide and 10 mm long and was segmented into 6 strips 
of 4 mm by 10 mm. The SSD covered a polar angle of 
= 90.0 ± 2.7° and an azimuthal angle of (f> = ±1.1° in 
the laboratory frame. The MCP assembly, consisting of 
two layers of MCPs and a single anode, had a sensitive 
area of 27 mm in diameter. Recoil carbon ions passed 
though the carbon foil and were detected with the SSD. 
The SSD provided the kinetic energy and the time-of- 
flight (TOF) from the target to the SSD. Secondary elec- 
trons, emitted from the carbon foil primarily by the pas- 
sage of recoil carbon ions, were accelerated and reflected 
to the MCP with electrostatic fields. The MCP detected 
those electrons and provided a better start time than the 
RF pulse from the AGS, but with efficiency too low to be 
used in the final analysis. Data collected with the MCPs 
were used in systematic checks of the measurement and 
a background estimation. 

Data were recorded for every AGS cycle. The primary 
trigger was a coincidence between the signal from one 
of the SSD strips and the RF pulse. The dead time of 
the data acquisition system was measured to be typically 
5 %. Energy calibration of each of the SSD strips was 
carried out with an 2 ||Am (5.486 MeV a particles) ra- 
dioactive source. 

Low energy recoil carbon ions were clearly seen as 
a kinematical correlation between E and Tssd, where 
E and Tssd are the energy deposit and the timing in 
the SSD strip respectively. Fig. || (a) shows the cor- 
relation obtained from one of the SSD strips. Energy 
losses in the target, the carbon foil and the SSD dead 
layer (22.6 ± 1.0 /ig/cm 2 Si equivalent) were corrected 
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to obtain the kinetic energy of recoil carbon ions. The 
detector resolution was determined from the time-energy 
correlation and the reconstructed mass distribution of re- 
coil carbon events. The resolution of Tggrj was 7.0±0.5 ns 
independent of E and consistent with the beam bunch 
width. The resolution of (XgsD — Tmcp), where Tmcp 
is the timing in the MCP, was 1.2 — 3.8 ns, decreasing 
with E. The energy resolution of the SSD, obtained from 
both the time resolutions and the mass distribution, was 
AE/E = 0.05/ ^JE (MeV) © C (quadratic sum), where 
the constant value of C varied in SSD strips from 0.05 to 
0.12. 

The momentum transfer in pC elastic scattering can 
be defined by the kinetic energy of the recoil carbon as 
tfi or by the TOF from the target to the SSD as tx- 
Taking into account detector resolution smearing, mo- 
mentum transfer was defined by the resolution-weighted 
mean, t = (wEtE + WTtT)/{wE + WT), where the weights, 
We and wt , were determined from the energy resolution 
and the time resolution respectively. The calculation was 
dominated by tE since we /{we + wt) ranged from 0.6 
to 0.9 depending upon t. The Grange covered to obtain 
An was 9.0 x 10~ 3 < -t < 4.1 x 10~ 2 (GeV/c) 2 . 

Recoil carbon events were identified by the mass re- 
constructed from the kinetic energy and the TOF from 
the target to the SSD. Fig. | (b) shows the mass distri- 
bution in the measured i-range. The mass resolution was 
2.1 — 3.8 GeV/c 2 as a function of t. A 2.5-sigma mass 
resolution cut was applied to select recoil carbon events. 
Backgrounds in the selected recoil carbon events were es- 
timated as a function of t. Backgrounds estimated from 
the target-empty data were 0.3 — 2.7 %. Backgrounds 
due to target fragments (dominated by a's) were esti- 
mated to be 0.8 — 11.9 % from the mass distribution 
reconstructed from the kinetic energy and the TOF from 
the carbon foil to the SSD, using (Tssd — Tmcp) with the 
mass resolution of 1.2 — 2.5 GeV/c 2 . Total background 
was 1.4— 13.9 % and was subtracted from the selected re- 
coil carbon events for each spin state separately. Angular 
distribution of the selected recoil carbon events for each 
t was obtained from the distribution of the hit position 
in the SSD. The recoil angle, 8, and the multiple scatter- 
ing width, A9, in the measured Grange were determined 
to be 89.0° < 6 < 89.7° and 0.8° < A0 < 1.3°. The 
recoil angle as a function of t was consistent with that of 
pC elastic scattering kinematics. A total of 2.2 x 10 7 pG 
elastic scattering events were thus identified for the An 
calculation. 

The analyzing power of pC elastic scattering was de- 
termined as a function of t using the formula , 

1 y/Nl-N^-y/Nj-N^ 
An = ■= , , (1) 

FB V^i • N R + V n l ■ N l 
where N\(N^) represent the number of events, in which 
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FIG. 3: The analyzing power, An, for pC elastic scattering in 
CNI region. The error bars on the data points are statistical 
only. The solid line is the fitted function from theory [[f| . The 
dotted lines are the 1-sigma error band of the fitting result. 
The dashed line is the theoretical function with no hadronic 
spin-flip amplitude (rs = 0). 

the scattered proton was in the left (right) side and the 
recoil carbon was in the right (left) side to the beam axis 
with the polarization direction up (down). Eq. (Q) can- 
cels out the acceptance difference between the left and 
right detector arm and the beam intensity difference be- 
tween the up and down polarized state. Each square root 
represents a geometric mean of two independent mea- 
surements of the same quantity which follows the rota- 
tional symmetry of the cross section around the beam 
axis. The resulting An is shown in Table | and Fig. || 
and is nonzero at the lower i-range with high statistical 
precision. The systematic error in the raw asymmetry 
measurement was estimated from a number of possible 
systematic effects, mainly from the uncertainty in the 
event selection of recoil carbon ions. The time-dependent 
variation of the measured raw asymmetry was consistent 
with a statistical fluctuation, x 2 1 dof =1.16. The normal- 
ization error in the Vb measurement was 12 %, which 
comes from the uncertainty of A p ^ . Our result of A n 
is in marked disagreement with the model-independent 
calculation with no hadronic spin-flip amplitude (rs = 0) 
shown in Fig. ||, 

The ratio of hadronic spin-flip to non-flip amplitude, 
rs, was determined by fitting the measured An with the 
theoretical formula of Ref. 0). Errors used for An in the 
fitting were the linear sum of three errors in Table, ffl. The 
fitted function with 1-sigma error band is shown in Fig. |^. 
The best fit gave a value of r$ as Rer$ = 0.088 ± 0.058 
and Imr 5 = -0.161 ± 0.226 with x 2 /dof = 0.57. The 
result of r5 and its associated x 2 contours are shown in 
Fig. | 

The magnitude of the hadronic spin-flip amplitude has 
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TABLE I: An as a function of t (GeV/c) 2 . Errors in t are the 
systematic error from the estimation of the energy loss at the 
target. The three errors in An are the statistical error, the 
systematic error in the raw asymmetry and that in the beam 
polarization respectively. 



-t ((GeV/c) 2 ) 
(0.90 ± 0.03) xlO" 2 
(1.26 ± 0.04) xlO -2 
(1.72 ± 0.04) xl0~ 2 
(2.20 ± 0.05) xlO -2 
(2.69 ± 0.05) xlO -2 
(3.17 ± 0.06) xlO" 2 
(3.65 ± 0.06) xlO' 2 
(4.13 ± 0.06) xlO -2 



At 



(1.30 ± 0.22 ± 0.35 ± 0.16) xlO -2 
(1.48 ± 0.11 ± 0.12 ± 0.18) xlO -2 
(1.19 ± 0.11 ± 0.17 ± 0.14)xl0~ 2 
(1.05 ± 0.13 ± 0.11 ± 0.13)xl0~ 2 
(0.12 ± 0.15 ± 0.18 ± 0.01) xlO" 2 
(-0.03 ± 0.18 ± 0.10 ± 0.00) xl0~ 2 
(-0.20 ± 0.21 ± 0.25 ± 0.02) xlO" 2 
(-0.15 ± 0.24 ± 0.29 ± 0.02) xlO" 2 




Rer 5 

FIG. 4: Ratio of hadronic spin-flip to non-flip amplitude, rs, 
and the 1-, 2- and 3-sigma contours of \ 2 ■ 



been a long-standing question. The E581/E704 experi- 
ment at Fermilab measured An for pp elastic scattering 
in CNI region at 200 GeV/c |]|. The result was con- 
sistent with no hadronic spin-flip amplitude within the 
limited statistical precision. Theoretical estimates of r§ 
based on several approaches have been made Q , indicat- 
ing |r 5 | < 0.15 from the E581/E704 result and \r 5 \ < 0.10 
for the RHIC energy range of ^ = 50 - 500 GeV. The 
present results using the carbon target provide a bet- 
ter determination of spin properties of scattering ampli- 
tudes which may help to reveal the dynamical mecha- 
nisms for scattering in the high energy asymptotic re- 
gion. To understand the asymptotic spin properties of 
scattering amplitudes more systematically, it is neces- 
sary to measure An in the CNI region at higher energies 
where the Pomeron amplitude is dominant. This may be 
performed by our method at RHIC. 

In conclusion, we measured for the first time the An 
for pC elastic scattering in CNI region of 9.0 x 10~ 3 < 



-t < 4.1 x 10~ 2 (GeV/c) 2 at 21.7 GeV/c with high sta- 
tistical precision. A nonzero value of r§ was obtained 
from A N to be Rer 5 = 0.088 ± 0.058 and Imr 5 = 
—0.161 ± 0.226. The experimental method used is ap- 
plicable to high energy proton polarimetry at RHIC Q . 
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